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Model sets (also called cut and project sets) are generalizations of lattices, and 
multi-component model sets are generalizations of lattices with colourings. In this 
paper, we study self-similarities of multi-component model sets. The main point 
may be simply summarized: whenever there is a self-similarity, there ars-ialso 
naturally related density functions. As in the case of ordinary model sets 0, we 
show that invariant densities exist and that they produce absolutely continuous 
invariant measures in internal space, these features now appearing in matrix form. 
We establish a close connection between the theory of invariant densities and the 
spectral theory of matrix continuous refinement operators. 



1 Introduction 

Mathematical quasicrystals are tilings or Delone sets obtained from the cut 
and project set or any other construction that ensures a pure point diffraction 
spectrum. In many examples, they are one-component sets, in the sense that 
their vertices (respectively their points) form just one translation class with 
respect to the limit translation moduleta - one well-known exception being the 
rhombic Penrose tiling whose vertices fall into four different classes. 

In view of the fact that practically all known physical quasicrystals are 
multi-atomic alloys, a systematic treatment of multi-component Delone sets is 
necessary and should give extra insight into the phenomena possible for non- 
periodically ordered structures, such as symmetries, inflation structure, etc. ICp 
make some progress here, one can use a generalized set-up of internal spaceE3, 
and/or study the mutual dependence of the components under inflations. 

Here, it is our aim to extend recent work on the existence of inflation in- 
variant measures on such sets from the one-component to the multi-component 
case. After introducing multi-component model sets and the notion of affine 
self-similarities for them, we show that the set of all affine self-similarities, 
based on a given linear inflation similarity, leads directly to the concept of an 
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invariant density on the points of the model set. Determining this invariant 
density necessitates moving the picture to the internal side where it becomes a 
problem of finding fixed points of a matrix refinement operator. We illustrate 
the nature of the solutions by examining the situation in a pair of examples 
based on the vertex set of a rhombic Penrose tiling, a 4-component model set. 

2 Multi-component model sets 

We begin with the notion of a normal cut-and-project scheme. By definition, 
this consists of a collection of spaces and mappings: 

R m R m x R n ^ R" 

U (1) 
L 

where R m and R n are two real spaces, n 1 and 7r 2 are the projection mappings 
onto them, and L C R m x R™ is a lattice. We assume that - is injective and 
that 7r 2 (I/) is dense in R n . We call R m (resp. R") the physical (resp. internal) 
space. We will assume that R m and R n are equipped with Euclidean metrics 
and that R™ x R™ is the orthogonal sum of the two spaces. For x lying in any 
of these spaces, |a;| denotes its Euclidean length. 
Let L := tt 1 (L) and let 

(.)* : L — > R™ (2) 

be the mapping ir 2 ° (t^i I ^ ) 1 - This mapping extends naturally to a mapping 
on the rational span QL of L, also denoted by (.)*. 

Let us now assume that L is a sublattice of another lattice M, also lying 
in R m+ ™. Necessarily, rank^(M) = m + n = rankz(Z), so [M : L] is finite. 
It is immediate that the pair (R m+n , M) gives rise to another cut and project 
scheme and a new Z-module M := ~k\ (M) of finite rank equal to m+n. Clearly, 
M lies in the rational span of L and [M : L) = [M : L] . 

We now choose r distinct | cosets iW := z l + L, r > 1, of L in M and r 
windows fl^\ . . . , Q\ r l inside the internal space R". Throughout this paper, 
we will make the following assumptions on all windows Q, in particular on the 
r windows that we have just chosen: 

Wl Q C R" is compact; 

'Distinctness is by no means necessary, and it might actually be useful for some applica- 
tions to drop this assumption which makes no difference to the later arguments. 
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W2 O = int(O); 

W3 The boundary of f2 has Lebesgue measure 0. 
With this data, we now define r subsets of R m by: 

A« := A^(fiW) := {x € | x* € flW } . (3) 

Since we are assuming that the cosets are distinct, the A^ are mutually 
disjoint. There are two ways in which we can view the collection of sets A^. 
The flat picture is 

r 

A := |JAW C M m (4) 

i=l 

in which we view the sets simultaneously in R m , as one would ultimately need 
for atomic models. The second, and mathematically more useful, way is the 
fibred picture in which we embed the sets AW into r copies of R m : 

A 1 := H A« c II R m , (5) 

i— 1 i— 1 

We call this fibred picture a multi- component model set. 

Let us pause to mention some of the nice properties of such point sets in 
the flat picture. 

Ml A is a Delone set, i.e. A is uniformly discrete and relatively dense. 

M2 A is a Meyer set: in addition to Ml there is a finite set F so that 
A — A C A — F. Equivalently, A — A is also a Delone set. 

M3 A has a well-defined density d, i.e. 

d = lim — — ^— - = lim — — — (6) 

vol(B s (0)) s^°o c m s m 

exists, where B s (0) := {x e R m | \x\ < s}, A s := A n B s (0), and 
c m '■= ^ s ^ ne volume of the unit sphere in R m . 

M4 A diffracts, i.e. it has a well-defined Bragg spectrum. 

These facts can be demonstrated by slightly rearranging the setup of 
the cut and project formalism and then appealing directly to the papers of 
Schlottmann 113 and Hof B where they are proved. Set G := M/L and let 
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a : M — > G be the canonical quotient map. We create a generalized cut and 
project scheme: 

W n < R m x ( R n x qj > R" X G 

U (7) 

N 

where N := {(i, a(ac)) | x <E M} c (M m x R") x G = R™ x (R™ x G) with 

7Ti(x) «- (a;,a(x)) -> (tt 2 (x), a(a;)). (8) 
Then, using the window 

r 

Q = |J (n^ x (9) 

J=l 

we recover A directly from (^). 

The limit in (||) is easily seen to be independent of the choice of origin for 
the Euclidean space; indeed, it even exists uniformly for sets: for any subset S 
of £1 with boundary of measure 0, the relative frequency of the points of (A s )* 
falling into S, as s — > oo, is vol(5)/vol(f2), and the convergence is uniform with 
respect to translations of the set S. See also Refs. 12, 13, and 11. 

Multicomponent model sets are a natural and somewhat parallel algebraic 
counterpart of tiling systems with various types of tiles or vertices. For ex- 
ample, in the case of the rhombic Penrose tilings (see example below), it is 
well known B that there are four classes of vertices corresponding to the four 
different coscts of the limit translation module (LTM) , the class of the LTM 
itself not being present. These four classes can be viewed as algebraic in origin, 
ultimately deriving from the structure of the prime ideal over 5 in the ring of 
cyclotomic integers for the fifth roots of unity, see also Ref. 1. 

As in the case of ordinary model sets, the generalization out of the domain 
of ordinary lattices entails the loss of translational symmetry. What often 
emerges instead, as one sees in many tiling models, is a very rich structure 
self-similarity. Our aim here is to begin to derive the analytical consequences 
of these new forms of symmetry. 

3 Self-similarities 

Definition 1 A self-similarity of the fibred model set A^ of raj is an r X r 

matrix t = tQ tV of maps 

t ji . a« — ► (10) 
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of the form 



t jl : x h-> Qx + v j \ 



(11) 



where Q is a fixed linear similarity map (i.e. a rotation followed by an inflation 
with a factor of q) stabilizing both L and M , and the v Jl lie in R m . 

We call Q the linear part of t and the constant q the inflation factor of t. 
Notice that this concept of self-similarity is not necessarily compatible with 
any sensible mapping, on M. m , of A into itself, unless the are independent 
of i,j. The latter situation happens, of course, in the one-component caseB. 

Note that Q naturally gives rise to an automorphism Q of the lattice M, 
i.e. an element of GLz(M), and a linear mapping Q* of K" that maps O into 
itself. From the arithmetic nature of Q we deduce that the eigenvalues of 
Q and Q* are algebraic integers and from the compactness of f2 that Q* is 
contractive. Furthermore, one can deduce that Q* is diagonalizable from 
the corresponding property of Q. Strictly speaking, we can only deduce that 
the eigenvalues of Q* do not exceed 1 in absolute value, but we will always 
assume that in fact they are less than one in absolute value. Such linear 
transformations induce so-called hyperbolic transformations on the associated 
torus M. m+n /M. Since Q also stabilizes L, Q stabilizes L and hence induces an 
automorphism of the group G = M / L. In the sequel, we will normally denote 
the contraction Q* by A to match other sources. 

To proceed, fix Q, once and for all, and look at the set T = Tq of all 
self-similarities t v = tQ jV of which have linear part Q. We define 

:= {tiel" Q*fiW + u c (12) 

We assume each ri JI is non-empty and satisfies the window conditions [W] 
above. (In fact, it is not crucial that the second window condition be fulfilled, 
as we will explain below, and we will use this flexibility in our examples. 
However, it makes the exposition a little less cluttered to assume it). 
We define 

L* := L+tf-Qz*) 

T Jl := {y G V 1 | y* G { ' 

It is straightforward to see that for any matrix r x r matrix of elements of M m 
we have 

t : = t v eT Q v ji G T ]i , for all j, i. (14) 

An invariant density on A^ is, by defintion, a set of non-negative functions 
pi : Jjii) _> 

R> , j = 1, • • • , r, satisfying the following properties: 
ID1 The function pi vanishes off A^^; 
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ID 2 The functions observe the equations 



pi(x) = | det(Q)| • lim £ £ p'W'H* ~ «)) i ( 15 ) 



ID3 The following limits exist: 



xeA 



where ^ = (^■ Ji ) is a non-zero matrix of non-negative numbers and w — 
(w , . . . , w r ) t 7^ is a vector (resp. an r x 1 matrix) of non-negative numbers. 
We will see that in order to solve for such invariant densities, it is necessary 
that the matrix v and the vector w satisfy: 

PF1 The spectral radius of v is 1, this is an eigenvalue of v, and w is the 
corresponding right Perron-Frobenius eigenvector. 

The matrix v encodes the degree of freedom to weight the points of one coset 
relative to another. In view of this interpretation, it would be reasonable to 
choose v to be a Markov matrix. In any case, the irreducibility of v, in the sense 
of non-negative matrices, would be sufficient (though by no means neccessary) 
to guarantee the additional property 

PF2 1 is a simple eigenvalue of v. 

The right eigenvector w (which has non-negative entries) may be normalized 
statistically, i.e. 53j=i w J = 1, whereupon it is uniquely determined. Note that 
such a normalization takes the relative densities of the points per unit volume 
into account, while one could also work with a normalization per point. 

In effect, the main condition ID2 says that, apart from a rescaling factor 
|det(Q)|, the density p->(x) at a point x of A^' is the averaged value of the 
densities of the points y which are mapped to x by an element of 7q, after 
these averages have been weighted by the non-negative matrix v. ID1 avoids 
the necessity of determing which of the possible preimages Q _1 (x — v) actually 
lie in in (|l5|), and ID3 is a normalization condition. 

4 Solving for the invariant densities 

We can solve for the invariant densities in the same manner as in Ref. 2. We 
assume the existence of a set of continuous functions p : flu'' — ► R> such 
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that supp(/ J ) C fi^ and f j (x*) = p 3 (x) for all x G A<- j \ We rewrite the 
equations [ID] on the window side. The uniform density of projection on the 
window side allows us to use Weyl's theoremllMJ to replace the averaged sums 
by integrals. To avoid undue proliferation of symbols, we revert again to using 
x as the variable, now living in R™. With A :— Q* , this leads to 

P(x) = |det(Q)| ^ =1 _^i_/^ ^(A-^x-z;))^; 

(17) 

w ° = /nco f j (x)dx. 

Setting Xs := V ^ S \ for any measurable subset S of R", we can rewrite the 
first equation in the form 

P(x) = |det(Q)| £V* / X nji (x-y)f {A^y)dy. (18) 

i— 1 ^ n 

Defining Y to be the matrix of functions (u^X^ji) and / = (f 1 , . . . , / r )*, we 
can write this more succinctly as 

/(a:) = |det(Q)| / Y(x-y)f(A- 1 y)dy (19) 

which expresses it in the form of a 'matrix convolution'. 

We can solve this system of integral equations by taking Fourier trans- 
forms, getting (with A* denoting the transpose of A) the matrix equation: 

/» = f(fc)./(A*fc). (20) 

This recursive formulation leads at once to the form of the desired solution (in 
Fourier space): 

oo 

/» = nw)')( fe )-/(°)- ( 2i ) 

£=0 

From k = in (pTj|), we get the consistency equation 

/(0) = w = f (0) • /(0) = i/ • /(0) , (22) 

so u; is an eigenvector for the matrix v = (v^ % ) with eigenvalue 1. This explains 
the necessity of the conditions [PF1]. 

The matrix linear operator R on the Banach space L^,(fl) of r x 1 column 
matrices of L 2 functions with support of the ith component in OW defined by 

R(F) = |det(Q)| / Y(x-y)F(A- 1 y)dy (23) 
7 



is an example of a matrix continuous refinement operator in the sense of Jiang 
and LeeB. We are looking for a 1-eigenfunction of this operator. The situation 
regarding the spectrum of R requires some care because there are two matrices 
involved: A and v. Under the assumptions [PF1] a solution g G C°°(f2) exists. 
If, in addition, [PF2] holds then there exists a unique C°° solution g (see Ref. 8, 
Theorem 3.1 and Remark 1 following it). In our particular case, it is not too 
hard to see that the infinite product (Elf) converges uniformly on compact sets 



to a vector of C°°-functions. Now, g also satisfies the defining equation (20) 
and g(0) is also a 1-eigenfunction of v, hence, after rescaling, can be taken to 
be w. Then it is straightforward to see that in fact g is our / of (pl|). Finally, 
taking the inverse Fourier transform of d2l]), we obtain 



g(x) = f(x) 



* (Y o A 1 ) 



(x)-w, (24) 



with obvious meaning to the convolution as applied to matrices. This is then 
the unique solution to the invariant density problem of ( |l7| ) . Defining 

p*(x) := f j (x*), j = l,...,r, (25) 

solves the invariant density problem on the multi-component model set A. The 
functions / J may be viewed as the Radon-Nikodym derivatives of a set of r 
absolutely continuous measures on the windows fiW . 



5 Example of Penrose tilings 

The simplest multicomponent model set of relevance to quasicrystal theory 
is the set of vertex points of the above mentioned phombic Penrose tiling. 
Recall the algebraic construction of such vertex sets El. Let £ := e 27 ™/ 5 and 
let Z[f] := {Ej=o m i£ j | r > 0,rrij G Z} C C. Clearly, Z[£] is the ring 
of integers of the cyclotomic field of the fifth roots of unity and is a rank 4 
Z-module with basis {1, £, £ 2 , £ 3 }. We can embed Z[£] into C 2 ~ R 4 by the 
mapping x t— > (x,x*), where * is the Galois automorphism Z[£] defined by 
sending £ to £ 2 . There is a homomorphism g : Z[£] — > Z/5Z defined by 
qQ2 m j£> J ) — 12 m j m od 5 , which is suitable to define cosets. 

Let P be the pentagon which is the convex hull of {1, £, . . . £ 4 } and define 
the windows ftW := P, := -P, il^ := tP, := ~tP. The notation 
is chosen to be compatible with the fact that g(r) — 3 mod 5. Define LW := 
{x G Z[£] | q(x) — i} and observe that . . . , are four different cosets 
of in Z[£]. We define 

:= {x e | x* G O^}. (26) 



8 



Then, A = [j A^' is the set of vertices of a (singular) Penrose tiling. More 
precisely, A is the union of 10 such singular sets, deviating from a regular set 
only in points of density 0, due to window condition [W3]. To obtain regular 
tilings, it is necessary to add a displacement 7 6 C to each of the windows 

i n (^6|), but ultimately this makes no difference to the measure theoretical 
considerations that are involved in the invariant densities on the windows, so 
we choose to suppress this additional complication here. 

We take, as our fixed linear self-similarity, the scaling maps Q = rlc, 
with corresponding contraction A = Q* = — -1q on the window side. The 
transition windows of @ become (with R 2 = C) 



{ueC 



(27) 



Determining these is a simple geometrical problem the outcome of which can 
be summarized in the following table 



(l/r 3 )P 
-P 

(l/r 3 + l/r)P 
-(l/r 2 )P 



{0} 
-(l/r 2 )P -(l/r)P 



(l/r)P 



(l/r 2 )P 
{0} 



(l/r 2 )P 
-(l/r 3 + l/r)P 
P 

-(l/r 3 )P 



(28) 



where P = QW. The ji entry of the table is the region fF\ We meet here 
a natural situation in which some of the windows are degenerate (the empty 
and singleton windows) and do not satisfy the window hypotheses. Since we 
rely on Weyl's theorem to connect the physical and internal sides of the pic- 
ture, and this is not applicable in such cases, we have to omit transition cases 
where vol(0 JI ) = 0. Otherwise, the normalized characteristic functions create 
"ghosts" , i.e. XQji is no longer a function but becomes a Dirac measure. We 
adopt the convention of choosing v Jl = whenever that happens. 

We are still free to choose the transition matrix v. For the purposes of 
illustration, we have chosen two examples and included figures that show the 
shapes of the corresponding invariant densities on the windows. 



5.1 Example 1 

In this example, we define 

/ vol(f^) \ 



ELivoi(^) 



/ (2-r)/4 (r-l)/4\ 

t/4 2-t t-1 (3-t)/4 

(3-r)/4 t-1 2-t t/4 

V(r-l)/4 (2-r)/4/ 



(29) 
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in which the contribution of points of type i to the density of points of type j 
is weighted by the areas of the corresponding transition windows. In this case, 
v is a Markov matrix (i.e. > and Y^j=i = 1)' so the PF eigenvalue is 
1, with corresponding eigenvector (0, 1, 1,0)*. We get non-vanishing densities 
only on the two windows fl^ and Q^ 3 \ a phenomenon caused by v being a 
reducible matrix (although condition [PF2] still holds). 




5.2 Example 2 
Here, wc define 
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V 2 








2/ 



(30) 
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This v is the transpose of a Markov matrix, whence the PF eigenvalue is still 
1, now with right eigenvector (1,1,1,1)*. This means that, on the 4 windows, 
the integrated densities must give equal values. 




-2 



Figure 2: Invariant densities on windows and QP^ 1 for Example 2. Notice that the two 
plots are scaled differently and that the invariant density on the right is defined on the larger 

of the two windows. 



In the Figures, we show the invariant densities on only two of the win- 
dows, a large and a small pentagon. Because of the inherent symmetry in our 
examples, the densities on the other two windows differ only in orientation. 
Though it is not necessary for v to be Markov, the use of a Markov matrix, 
which preserves probabilities, seems appropriate for potential applications of 
this picture to a probabilistic version with positive entropy densityCrlj. Furthc 
examples with different scaling factors, in particular with q = t? where the 
invariant densities flatten out at the top, will be given elsewhereS. 
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